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NONLINEAR PROBLEMS IN FLIGHT DYNAMICS INVOLVING AERODYNAMIC BIFURCATIONS 


Murray Tobak and Gary T. Chapman 
Research Scientists 

NASA Ames Research Center, Moffett Field, California 9*1035, U.S.A. 


SUMMARY 

Aerodynamic bifurcation is defined as the replacement of an unstable equilibrium flow by a new stable 
equilibrium flow at a critical value of a parameter, A mathematical model of the aerodynamic contribution 
to the aircraft’s equations of motion is amended to accommodate aerodynamic bifurcations. Important 
bifurcations such as, for example, the onset of large-scale vortex-shedding are defined. The amended 
mathematical model is capable of incorporating various forms of aerodynamic responses, including those 
associated with dynamic stall of airfoils. 


1 . INTRODUCTION 

In recent years we have become increasingly aware of the connections which link our approach to non- 
linear problems in flight dynamics to the body of theory underlying the rapidly growing field known as 
"nonlinear dynamical systems." In Ref. 1, we attempted to identify the connections explicitly within a 
framework designed to facilitate study of the aerodynamic contribution to the flight-dynamics system. The 
framework is composed of parallel observational and mathematical components. A study of observations of 
fluid flows singled out four important elements for the observational component. First, flow patterns 
have definite structure (e.g., three-dimensional vortical structures). Second, structures undergo sys- 
tematic changes with variations in parameters such as angle of attack or Reynolds number. Third, under 
some conditions the changes lead to chaos (e.g., turbulence). Fourth, the structures have definite spa- 
tial and temporal scales . Parallel elements of the mathematical component feature the language of topol- 
ogy and concepts drawn from dynamical-systems theory. The premises of the theoretical framework are: 

(1) Structures are describable in topological terms. (2) Changes in structure are describable by 
bifurcation theory. (3) Chaos is describable by the theory of strange attractors and fractals. (4) 

Scales are describable by group theory ideas. Recently (Ref. 2), we proposed that the framework may be 
useful as well for the understanding and the potential modeling of turbulent flows. 

Here, we shall focus on one element of the framework, namely, change, and its representation by means 
of bifurcation theory. Bifurcation theory has two principal realms of application in flight-dynamics 
studies. First, when the form of the aerodynamic contribution to the aircraft’s equations of motion is 
known, and the equations of motion are of the appropriate form, bifurcation theory provides a useful means 
of analyzing the aircraft’s motion, particularly in the vicinity of stability boundaries. This first 
realm was the subject of study in Ref. 3, and the study is pursued further in Ref. 4 at this conference. 
The second realm, which will be the subject of this paper, is concerned with the mathematical modeling of 
the aerodynamic contribution itself. 

Reference 5, a companion paper of Ref. 1, was a reconsideration of an approach to mathematical model- 
ing of the aerodynamic characteristics, in which nonlinear indicial responses are used as a basis for 
arriving at the form of the aerodynamic response to arbitrary motions. Without modification, validity of 
the approach rests on a fundamental assumption that the indicial response exists and is unique in every 
incremental step of the motion. The question, "When does the fundamental assumption fail?" led to the 
recognition that the assumption may be invalidated through the mechanism of flow instability. On exceed- 
ing a critical value of a parameter (e.g., angle of attack), the flow field corresponding to the equilib- 
rium state of an aerodynamic response may become unstable to small disturbances. The replacement of the 
unstable state by a new stable equilibrium state satisfying the same boundary conditions implies failure 
of the uniqueness condition, and may imply failure of the existence condition as well. Thus, the equilib- 
rium state undergoes a bifurcation, and the circumstances are precisely those that bifurcation theory is 
designed to address. Although the requisite amendments were not carried out in Ref. 5, it was concluded 
that bifurcation theory would provide the means to extend the mathematical model so as to acknowledge the 
occurrence of aerodynamic bifurcation phenomena. As a consequence, the model would have the potential of 
accounting for those critical points in maneuvers where sudden and dramatic changes in flow structure may 
occur. 

In this paper, we shall try to validate the conclusion of Ref. 5, and demonstrate how the mathemati- 
cal model can be amended to accommodate the occurrence of aerodynamic bifurcations. A number of important 
aerodynamic bifurcations will be defined and some of their potential effects on dynamical systems illus- 
trated. Following the analysis presented in Ref. 5, we shall show how the approach to the modeling 
involving nonlinear indicial responses can be made compatible with an approach based on the use of non- 
linear algebraic functional expansions. The latter approach, which has the virtue of concreteness, was 
formulated originally (Ref. 6) to address nonlinear problems in the field of electrical circuits and sys- 
tems. Its adaptation here exemplifies the community of interests developing across various disciplines 
through the unifying medium of research in nonlinear dynamical systems. 
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2. MODELING INFLUENCE OF PAST MOTION BY PULSES 

The following is a slightly modified version of the analysis originally presented in Ref. 5. Again, 
for clarity, the two-dimensional wing is adopted as illustration, although the results will have more 
general bearing. 

Let the wing move away from a coordinate system whose origin is fixed in space at the center of grav- 
ity at a time 5 - 0. The distance traveled by the center of gravity along the flightpath is measured by 
a coordinate s. Let the center of gravity move at constant velocity V Q , so that the trace of its path, 
plotted against time 5, i 3 a straight line. This is shown on Fig. 1 . The wing is allowed to undergo 
changes only in the angle of attack a, where a is the angle between the velocity vector and the wing 
chord line. Projections of the leading and trailing edges of the wing onto the plane containing the 
velocity vector are maximum when a - 0. These maximum projections also trace out straight lines on 
Fig. 1 , parallel to the trace of the center of gravity. 

As illustrated in Fig. 1 , let the angle of attack a be zero for all time 5 except at £ - ^ , 
where a pulse occurs of amplitude aC^) and of infinitesimal duration . Consider the response to 
the pulse at ^ of the lift coefficient C L at a measuring time t subsequent to ^ . The lift 
response at t will be a function of the elapsed time t - ^ and the amplitude of the pulse a(^ 1 ). If 
we assume that there will be a range of a (0 < a < ) in which the dependence of C L on aC^) is ana- 

lytic, then we can write the response in C L as a Taylor-series expansion in aCf^) of the form 

c L (t> - c L (t - ■ 2 a n (t " V [o( V ]n A5 i (1) 

n 

The first of the forms in Eq. (1) will be used in the subsequent analysis to distinguish between direct 
(subscript dir) and interference (subscript lnt) effects. 

Now let us consider the response in C L at the measuring time t to a pair of pulses located at 
and £ 2 with 5 ,5 2 < t. Here, in addition to the direct influence of each of the pulses acting as if in 
isolation, the interference between the pulses also will influence the lift. The interference effect can 
be written in a form resembling a product of responses to single pulses 

C L ■ 2 t> mn (t - e,.t - e 2 )[a(e i )]“’[a(C 2 )] n AC, A£ 2 (2) 

int ,2 m,n 

where the subscript (int, 2) means "interference between a pair of pulses." The form vanishes properly 
with the vanishing of either of the pulses and retains the analytic dependence on angle of attack. With 
the addition of the direct Influence of the two pulses, the lift coefficient at time t takes the form 


C L (t) - C L (t - C^aCC,))^ * C L (t - e 2 .aU 2 )) dir ♦ C L (t) lntf2 (3) 

The process of adding pulses can be continued indefinitely in the same way. At the next stage, the 
interference between triplets of pulse3 must be considered as well as that between pairs. Going to the 
limit of a continuous distribution of pulses starting at time 5-0 yields a summation of multiple inte- 
grals having the form 


C. (t) - C. + C. + C. 

L dir int, 2 L Int,3 


with 


C L - E I a (t - 5 1 )[a(5.)] n d5. 
dir n 

- 2 «2 r 2 b mn (t “ “ ^ 2 > ) ]" 

i nt , 2 m , n o o 

L,„. ■ E /' “V 1 " «3 jf 3 “V** «a 

lnt, 3 m,n,p o o 

/'* °mnp (t - ~ e 3 )Ca(C 1 )] P <*, 


w 

(5) 

( 6 ) 

(7) 


Written as a nonlinear functional expansion, Eq. (4) represents the lift coefficient at time t in 
response to an arbitrary variation of a over the time interval zero to t. The form of Eq. (4) confirms 
an important point made in Ref. 6. It will be noted that a partial summation consisting of the leading 
term from each of Eqs. (5), (6), . . . forms a Vol terra series (Refs. 6,7). The fact that there are terms 
left over confirms that the a priori adoption of a Volterra series to represent the lift coefficient would 
have been insufficiently general to accommodate the Taylor series form of the dependence on angle of 
attack . 
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3. FORMATION OF INDICIAL RESPONSE 


Given Eq. (4), one can now use it to form the indicial response in accordance with the definition in 
Ref. 8. To indicate the form of the result, it will suffice to consider terms in Eq. (4) only through the 

series representing . Two motions need to be considered. In the first, the wing undergoes the 

i n t , 2 

motion under study ot(£) from time zero up to a time £ - x, where x < t. Subsequent to t, a is held 

constant at a(x). Thus, in Eqs. (4), (5), and (6), 


a 1 (£) - a(0 ; 
- a(x) ; 


0 < £ < T 

e * t 


( 8 ) 


The direct and interference contributions to C L (t) take the form 

C ■ E f 1 a (t - e 1 )[aU 1 )] n d£, ♦ £ C«(T)] n C a (t - 5. 
dir n * f o n *x 1 

■5. 


)dC, 




“ L ^ I .o.(C 2 )] m d5 2 f 2 b (t - 5,.t - 5,)[a(e,)] n <35, 

int , 2 ra,n "o o 


*r j i 


E / T ta( 

n,n ¥ o 

* E [a(T)] m C d? 2 C b (t - C,.t - e,)[ aU )] n d? 

m , n *0 

* E ca(t) m+n ] f d^ 2 / ?2 b mn (t - c,.t - e 2 )d 5l 


(9) 


m,n w t "x 

In functional notation (Ref. 7), the lift response to the same motion is 


C L (t) - C L [ ai U)] - C L Ca(0;t,i] (10) 

where the notation in the first form should be understood to mean that C L (t) is a functional of the vari- 
ation 0^(0 over the interval 0 < £ < t. The second form, which is defined in Ref. 8, is designed to 
make more explicit that represents the motion under study a(0 over the interval 0 < £ < x, but 

is constrained to remain constant at a(x) for x £ £ £ t. Thus, C^tt) is both a functional of a(£) and 
a function of t and x. 

In the second motion, the wing undergoes the same angle-of -attack history a(£) up to time x. Sub- 
sequent to x, the angle of attack Is again held constant, but is given an incremental step change Act 
over its previous value of a(x). Thus, in the second motion. 


a 2 U) - a(0 ; 

- a(x) ♦ Aot ; 


0 < K < x 

e * t 


} 


(id 


The direct and interference contributions to C L (t) become 

c, - E f T a n <t - e.KaU )] n d? ♦ E Mt) + Aa] n C a (t - e )d£ 

u dir n *0 " 1 


^int^ m,n ¥ o 


E f Co(£ 2 )] m d? 2 f ^ 2 b (t - e..t - d5, 

m,n •'o *o 

+ E £“<*) + A “] m / d ^ 2 f *> (t - 5 r t - 5 2 )Ca(£ 1 )] n d£. 

m,n *x * 

+ E + Aa] m+n j d5 2 J ^ 2 b^tt - 5,,t - 5 2 )de, 


ra,n 

Again, written as a functional, the lift response to the motion a 2 (£) is 

c L (t) - C L [a 2 U)] 


( 12 ) 


(13) 


The indicial lift response is formed by taking the difference between Eqs. (9) and (12) (or between 
Eqs. (10) and (13)), dividing by the incremental step Aa, and going to the limit as Act ■* 0. Existence 
of the limit is ensured by the analyticity of the functional dependence on a(£). The result is 
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lira 

Aa-K) 


AC L (t) 

Aa 


lira 

Aa-K) 


fC L ta 2 U)l “ C L [o 1 ( 03 ^ 
l Aa 


- £n[a(T)] n 1 
n 


s: 


a n ( t - 


♦ E "'[“(t)]" 1 ' 1 J d e 2 b mn (t - 5,.t - C 2 )Ca(C,)] n d5, 


+ E (m + n)[a(x)] m+n ' 1 f d? f^ 2 b (t - (.,t - !• )d£, (lit) 

m,n x t 

Eq. (14) reveals the form of the indicial lift response to a step change in angle of attack in terras of 
functional expansions. The first and third terms on the right-hand side of Eq. (14) do not depend on the 
past motion, but only on the level of the angle of attack at which the 3tep was made. The second term 
depends on the past motion, since a(^ 1 ) with 0 < C < t appears within the integral. The leading term 
of this past dependence has the form 


/ ^ r b n (t - 5 r t - 5 2 )a(5,)d5, 


Dependence on the past thus arises from interference effects between pulses prior to t» the origin of the 
step, and perturbation pulses of O(Aa) originating subsequent to x. In the general case, and just as 
before (Ref. 8), the indicial response is itself a functional. 


In fact, in formal terms (see Ref. 9 for a useful account of the essentials) the operations involved 
in forming the indicial response amount to taking a Fr£chet derivative of the functional C^Ca^S)]. To 
see this, we adopt the notation of Ref. 9, letting 


h 


en 


where we set 


so that 


e - Aa 
n - 0 ; 
- 1 ; 


o < e < x 
e * t 


a 2 U) - a } (0 + en 


Then, following Ref. 9, we have 


lira 

Aa-K) 


AC L (t) 

Aa 


lira 

e+0 


{£) ♦ en] - 


(15) 


(16) 


(17) 


-iF c L Ca i (t) + en] l e -o 

-c^[a,(5)]n (18) 

The equivalence of the operations yielding Eqs. (14) and (18) should be evident. Finally, to conform with 
the notation of Ref. 8, we set 

AC (t) 

lira — C. , [a 1 (5)]n ■ C. [a(5);t,x] (19) 

. _ fla L» i L 

Aa-K) a 

where the indication of a separate dependence on t and x rather than on elapsed time t - x alone 

should be noted. It can be easily verified that the first and third terms in Eq. (14), which depend only 

on a(x), indicate a dependence on t - x alone; however, as a consequence of its dependence on the past 

motion, the second term cannot be ca3t as a function of t - x alone. 


4. GENERALIZED SUPERPOSITION INTEGRAL 


in 


Just as before, Eq. (19) can be used to form a generalized superposition integral for the response 
C L to an arbitrary angle-of -attack variation. The result is 


c L (t) 


C L (t;a(0)) 


f 


C L [a(5)it,x] — dx 


(20) 
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with 


C L (tsa(0)) ■ C L [a(£);t,0] (21) 

It is important to understand the meaning of the first term in Eq. (20). It is the contribution to the 
lift at time t caused by the flow field existent at the initial instant £ - 0. It would be the only 
contribution if the angle-of-attack variation subsequent to £ - 0 had remained constant at a(0). A 
reading of the notation of the second form for this term in Eq. (21) (cf. the definitions following 
Eq. (10)) provides a more complete interpretation: The first term in Eq. (20) is the lift at time t due 

to an angle-of-attack variation a(£) which is equal to the motion history prior to 5-0, and which is 
constrained at £ = x - 0 to remain constant at a(0) for all £ £ x - 0. by substituting Eq. (1*1) for 
[a(£) ; t , x] in the integral term in Eq. (20) and carrying out the integration, one will verify that the 

form of Eq. (H) is restored through terms of the series representing C T (i.e., through Eqs. (5) 

int,2 

and (6)). We conclude first that an approach to modeling based on nonlinear indicial responses and gener- 
alized superposition integrals is compatible with one based on nonlinear functional expansions. Second, 
validity of the approach rests on the assumption of FrSchet differentiability of the lift response 
C L [ot(£) ;t, t] over the interval 0 < x < t. Invalidation of this assumption may signal the occurrence of 
an aerodynamic bifurcation. 


5. CONNECTION WITH BIFURCATION THEORY 


Consider the first of the two maneuvers involved in the formation of the indicial response, that is, 
the maneuver (cf. Eqs. (8) and (10)) to which the lift response C^EcKO ;t,x] corresponds. The angle of 
attack attains a constant value ot(x) subsequent to t, and it is reasonable to expect that the flow field 

at the subsequent time t will approach an equilibrium state that corresponds to this fixed boundary 

condition as the elapsed time t - t •* ». Heretofore (Ref. 8), the analysis ha3 contained an implicit 
assumption that a3 the flow field approached the equilibrium state It became time-invariant, which meant 
that the corresponding lift coefficient C L approached a unique constant value C L (~;a(x)) as 

t - t -*> ®. As long as this was true, it was reasonable to expect that an incremental change in a(x) of 

0( Aa) would result in an incremental change in C L (*;ct(x)), likewise of 0(Aa). We now recognize that at 
least two conditions must be satisfied in order that this be true: (1) The change of 0(Aa) must not 
result in a loss of analyticlty in the dependence of C L [a(£);t,x] on a(£). (2) The time-invariant equi- 
librium flow represented by C L («;a(x)) must be asymptotically stable to small perturbations. We deal 
with the second condition here, assuming for the moment that the first condition has been satisfied. 
Although exceptions to sufficiency can be envisioned, we shall restrict our study to the class of flows 
for which asymptotic stability of the equilibrium state i3 not only necessary, but also sufficient, to 
ensure that 0 L [a(O;t,x] will be Fr^chet differentiable. 

It can happen, however, that as a(x) is increased in small increments, a critical value of a(x) can 
be reached at which the stationary equilibrium flow represented by C L (°°;a(x)) will no longer be stable to 
small perturbations. Loss of asymptotic stability signals loss of Frechet differentiability at the criti- 
cal value of a(x). With the stationary equilibrium state asymptotically unstable, the system will seek a 
new equilibrium state, and of those that are possible will select one that can remain stable to small 
perturbations. This is precisely the situation that bifurcation theory Is designed to address. Bifurca- 
tion theory provides a means of classifying and characterizing the properties of the new equilibrium 
states that can arise when the given equilibrium state becomes unstable. We call the replacement of an 
unstable equilibrium state by a new stable equilibrium state an aerodynamic bifurcation. We will now show 
how the mathematical model Eq. (20) can be amended to acknowledge its occurrence. 


Referring to Fig. 2, let us assume that the angle-of-attack variation under study a(x) passes 
through a critical value at a value of x - x q within the Interval of interest 0 < x < t. If the 

angle-of-attack variation a(x) were constrained to remain constant at a Q for all x > x Q , the corre- 
sponding lift response C^EatS) ;t , x ] would begin to depart from an initially infinitesimally close neigh- 
boring response, such that C L Ea(£) ;t , x ] would not be Frechet differentiable. Thus, we cannot allow the 
integration in Eq. (20) to proceed beyond x^ without acknowledging the loss of Frechet 
differentiability. Let us assume that for all other angles of attack within the range a(0) < a < a(t) 
the lift responses represent asymptotically stable flows permitting Frechet differentiation. Then we can 
isolate the critical state by stopping the integration in Eq. (20) just short of x Q and starting again 
just beyond . Within the isolated interval - e < t < + e, we acknowledge the aerodynamic 

bifurcation by allowing the solution to change discretely to a new equilibrium state. Thus, 

c L (t) - C L (t;a(0)) ♦ f° C L [a(E) it , t] ^ dT 
J o a 


where 


f C L [a(0;t,T] ^7 dT + AC L^ t;a c^ 

"x +e a 
c 


( 22 ) 


AC L (t;a c ) “ c L E a <S);t,x c + e] - C L Co(5);t.T c - e] (23) 

This procedure was introduced earlier (Ref. 8) to acknowledge the occurrence of mean or constant 
values of the lift coefficient in the equilibrium state C L (“;a(x)) which could exhibit the discontinuous 
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and double-valued behavior characteristic of hysteresis. This was a case in which the first of the two 
conditions just stipulated, namely, analytic dependence of C L [a(£);t,t] on a, was violated. The partic- 
ular way in which this occurred is common in aerodynamics: The variation with a of the lift coefficient 

in the equilibrium state C^(«;a(t)) develops a fold at a critical value of a, a c , so that the slope of 
C L («;a(t)) versus a becomes infinite there. A jump in C^(®;a(i)) necessarily ensues to the upper 
branch of the folded curve with an infinitesimal increase in a beyond a , and hysteresis follows on the 
return route. We show this as sketch (a) in Fig. 3 along with, for comparison, sketches (b) and (c) of 
the two most common forms of aerodynamic bifurcation. These sketches represent bifurcation to another 
time-invariant equilibrium state. The distinction between the class represented by sketch (a) (the class 
studied in Ref. 8) and the examples of aerodynamic bifurcations should be noted. In sketch (a) there is 
only one solution curve; no new branches resulting from instability appear. Discontinuous and double- 
valued behavior result merely from the fold that develops in the solution curve. In the case of super- 
critical bifurcation, sketch (b), a new branch of stable solutions appears which is connected to the old 

branch at the critical point a Q . Beyond a c , two branches of solutions exist satisfying the same bound- 

ary conditions, of which one is stable and the other Is unstable. Supercritical bifurcation solutions 
start from zero and are initially infinitesimal as a increases beyond ct c , so that Eq. (23) would be 
zero. Supercritical bifurcation is not an important mechanism by which time-invariant equilibrium states 
exchange stability. The bifurcation can be easily superseded by smooth alternative solution curves 
through a variety of means (cf. Ref3. 10-12). In subcritical bifurcation, sketch (c), the bifurcation 
solution curve that branches away from the now unstable branch is itself unstable. If it turns back, as 

indicated in sketch (c), the upper branch will be stable. A discontinuous Jump in C L (®;a(x)) results 

a3 a increases beyond a c » and hysteresis follows on the return route, just as in sketch (a). Here, of 
course, Eq. (23) will have a value. 

Figure 3 has had the important function of reminding us that there are at least two mechanisms which 
can invalidate the conditions permitting Fr^chet differentiation of the lift response C^tciU) ;t, t] . 

There may be others. Joseph (Refs. 11,12) notes, for example, that isolated solution branches of flows in 
the equilibrium state can exist which are not connected to other solution branches through bifurcation. 
Later, we shall entertain the possibility that the initial flow as well as the stability of flow in the 
equilibrium state may influence Fr^chet differentiability. A theory of the means of invalidating Fr£chet 
differentiability of the lift response offers the possibility of an alternative to bifurcation theory. 
Inasmuch as it would incorporate the latter. It would be a potentially more inclusive way of classifying 
aerodynamic phenomena that are important In flight-dynamics applications. 


6. QUALIFICATIONS 

Before proceeding with the study of aerodynamic bifurcations, we must introduce some important quali- 
fications. Our current understanding of fluid-dynamic bifurcations stems almost entirely from results of 
analytical and computational studies of laminar flows governed by the incompressible Navier-Stokes equa- 
tions. Within these confines, studies have been further limited to flows in bounded domains. Boundary 
conditions have been simple enough to permit reduction of the linear stability problem (the necessary 
first step in the study of bifurcations) to the form of ordinary differential equations (e.g., the Orr- 
Sommerfeld equations). The principal findings on which our current understanding is based are summarized 
in Refs. 10-13. Qualifications are necessary when we attempt to extrapolate this body of knowledge to the 
study of the aerodynamic bifurcations that occur under full-scale flight conditions. These involve turbu- 
lent compressible flows which, if averaging and modeling are Introduced, will be governed not by the com- 
pressible Navier-Stokes equations but by evolution equations of unknown forms. Further, the domains of 
the flows are typically unbounded , rather than bounded. Finally, the boundary conditions, reflecting the 
geometries typical of modern aircraft, ordinarily will not allow reduction of the linear stability problem 
to the level of ordinary differential equations. 

The situation requires us to stipulate a list of presuppositions which will be the basis for the 
extrapolations we make. The list is inspired by one that Lanford compiled to justify the study of turbu- 
lence within a framework based on the incompressible Navier-Stokes equations (Ref. 14). Our list will 
merely supplement Lanford' s, taking additional account of turbulence modeling, compressibility, and the 
necessity to consider flows in unbounded domains. The issue of boundary conditions is not addressed on 
the assumption that their increased complexity does no more than increase the role of numerical 
computations. 

6.1 Turbulence Modeling 

Let us consider the flow field to which a typical lift response C L [a(5) ; t , x] corresponds. We assume 
that the flow field at time £ - x is given (to within a certain accuracy), and we must determine the 
subsequent flow at time t. Inasmuch as a(£) remains constant at ct(x) for all 5 - t > x, we have fixed 
boundary conditions. This form of the fluid-dynamics problem 13 known as an Initial boundary value prob- 
lem ( IBVP ; cf. Ref. 11). We assume that any given realization of the IBVP is governed with sufficient 
accuracy by the compressible Navier-Stokes equations. Paralleling Ref. 14, in the absence of proof we 
simply postulate that a solution of the IBVP would exist for all £ - t fc x and would be unique. How- 
ever, the Reynolds number is sufficiently high so that a solution would exhibit chaotic behavior. Implying 
a sensitive dependence on the Initial conditions. Thus, small perturbations of the initial conditions 
would cause repetitions of the given motion to result in a set of motions having stochastic properties. 
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The situation calls for the application of some form of averaging of the compressible Navier-Stokes 
equations to suppress the appearance of the small-scale chaotic structures. After averaging, the reper- 
tory of solutions of the resulting equations rau3t remain faithful to that of the unaveraged equations, 
neither allowing uncharacteristic solution forms nor losing essential ones. The necessity of closing the 
averaged equations in conformity with this requirement introduces the turbulence modeling problem. The 
convergence of interests this situation represents for the fields of research concerned with turbulence 
modeling and with the modeling of aerodynamic responses ha3 been noted in Ref. 5. In Ref. 2, we have 

tried to show how the theoretical framework underlying the deterministic approach to the study of turbu- 

lence might assist the modeling effort, particularly at this juncture, where it is important to give ade- 
quate representation of bifurcation phenomena. We cannot go further here in specifying modeling require- 
ments except to postulate the one feature that must be retained in common with the original Navier-Stokes 
equations: Whatever the form of modeling adopted, it must remain possible to cast the modeled equations 

in the form of autonomous evolution equations, namely, 

u t - HCu.JD (21)) 

where u is the velocity vector and X is a parameter (e.g., angle of attack, Reynolds number, Mach 
number). We believe that the premise is a reasonable one. Involving no more in principle than the use of 
orthogonal projections (Refs. 15,16) to eliminate gradient terms (e.g., the pressure term) from the mod- 
eled equations. 

6.2 Compressibility 

We shall assume that modeled equations of turbulent flows resulting from averaging the compressible 
Navier-Stokes equations remain at least as well-behaved as the original equations and their incompressible 

counterparts. Of course, the new parameters introduced by compressibility (e.g., Mach number) will them- 

selves have critical values at which additional bifurcations should be expected. 

6.3 Unbounded Domains 

Let us return to the IBVP represented by C L [a(£);t,T] and assume that the premises of the previous 
sections are now in force. Henceforth, the equations governing the flow field are assumed to be modeled 
equations of turbulent flow, of known form, and satisfying the form of Eq. (24) by suitable projections. 

The issue concerning the absence of bounds on the flow domain arises in evaluating the stability of 
the flow field in the equilibrium state. We say that the equilibrium state is reached In the limit as 

t - t + ", so that transient effects associated with the initial conditions have vanished. A study of the 

stability of the equilibrium flow must be undertaken first, since bifurcation of the equilibrium state 
will not occur unless a critical value of the parameter in question (here, a) exists at which the flow 
field in the equilibrium state becomes unstable. Stability is evaluated by determining whether an infini- 
tesimal perturbation of the equilibrium flow decays or grows with time. The perturbation’s fate can be 
determined from a linear spectral problem, the governing equations being obtained from a linearized per- 
turbation of the equations governing the equilibrium flow. Typically, with a bounded domain, It can be 
shown that solution of the spectral problem yields a countably infinite set of isolated eigenvalues in a 

complex plane. Stability is indicated if all of the eigenvalues lie in the left-hand plane. On the other 

hand, little is known about the spectrum of eigenvalues for unbounded domains. In the case of flows in 
unbounded domains governed by the incompressible Navier-Stokes equations, it is believed (Refs. 17,18) 
that an incomplete set of discrete eigenvalues exists which again controls stability, and that there is 
additionally a continuous spectrum having negative real part (i.e., lying in the left half of the complex 
plane). We postulate that what is believed to be the case for flows in unbounded domains governed by the 
incompressible Navier-Stokes equations is in fact the case for our problem. 

6.4 Conclusions 

Accepting the last premise (which implies accepting all of the previous ones) leads to several con- 
clusions concerning the role of bifurcation theory in the modeling of aerodynamic responses. First, since 
we have postulated that stability of the equilibrium state continues to be controlled by the discrete part 
of the eigenvalue spectrum alone, entailing no change from the rules governing bounded domains, we con- 
clude that the role of bifurcation theory regarding the equilibrium state should carry over intact to - 
flows in unbounded domains. That is, all of the Instability mechanisms that have been uncovered by stud- 
ies of flows in bounded domains governed by the incompressible Navier-Stokes equations should have coun- 
terparts in flows in unbounded domains governed by modeled turbulent flow equations. With the onset of 
instability of the equilibrium state, bifurcation theory again should be capable of classifying and char- 
acterizing the properties of the new stable equilibrium states that can arise to replace the unstable one. 

On the other hand, concerned as it is only with the equilibrium state, bifurcation theory does not 
suffice to completely resolve the IBVP from which the indicial lift response is derived. We note that 
when the eigenvalue spectrum lies in the left-half plane, indicating stability of the equilibrium state, 
the IBVP for an infinitesimal step perturbation of the boundary condition a(£); £ £ t is a linear 
problem. In fact, the transient flow field from which the indicial response [a(5);t,x] is derived can 

be obtained from a suitable superposition of the eigensolutions (including both She discrete and continu- 
ous parts) of the linear spectral problem. The superposition is made determinate by the necessity of 
matching a prescribed flow field at the initial Instant £ - t. This Is how the dependence on the pa3t 
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motion is expressed, which makes C L [a(£);t,T] a nonlinear functional even though derivable from linear- 
ized perturbation equations. Thu3, recalling the derivation by means of functional expansions, we affirm 
that so long as stability of the equilibrium state ensures that C L [a(£);t # T] will be Frdchet differenti- 
able, the differential AaC^ [a(£);t,*r] can be determined from a linear problem. With the onset of insta- 
bility, however, the perturbation flow no longer can be obtained from the solution of the linear spectral 
problem, since at least one of the discrete eigensolutions would grow without bounds. Thus, the incremen- 
tal change in lift coefficient given by Eq. (23) must be determined from the fully nonlinear perturbation 
flow equations. 


7. AERODYNAMIC BIFURCATIONS 

In Ref. 1, we postulated six major subdivisions in the form of aerodynamic force and moment responses 
by means of a set of sketches, and these are reproduced in Fig. 4. Although the ordinate f a in each of 
the sketches can be made to represent a variety of characteristics, it will be convenient here to let 
f a represent the lift coefficient in the equilibrium state. We now briefly outline how modeling of the 
aerodynamic responses incorporating the six forms of lift coefficient may be effected by means of the 
amendment we have introduced to accommodate aerodynamic bifurcations. 

7.1 Major Subdivisions 

First, it should be clear that aerodynamic systems represented by sketches (a) and (b) in Fig. 4, 
involving linear or nonlinear single-valued lift variations, do not require modification of the treatment 
contained in, for example. Ref. 8. In formal terms, these systems can be characterized as having aerody- 
namic responses C L [a(£);t,i] that are Fr^chet differentiable over the entire angle-of-attack range of 
interest . 

Sketch (c) in Fig. 4 introduces a class of multivalued solutions of the flow in the equilibrium 
state. Mean or constant values of the lift coefficient in the equilibrium state exhibit the discontinuous 
and double-valued behavior characteristic of hysteresis. As we have already noted, such cases were 
treated in Ref. 8 by a procedure similar to the one described here. The class includes both solution 
curves having folds, and subcritical exchanges of stability between time-invariant equilibrium states. 

Sketches (d) and (e) in Fig. 4 introduce bifurcations that result in time-dependent equilibrium 
states. Accommodating these phenomena within the mathematical model will require the use of Eqs. (22) 
and (23). Hopf bifurcation is indicated in sketch (d), wherein a formerly stable time-invariant equilib- 
rium state is replaced by a time-varying periodic equilibrium state. The difference in lift responses 
resulting from these two forms of equilibrium states is reflected In Eq. (23). Hopf bifurcations are 
common in aerodynamics. Physically, the usual origin of such a large-scale oscillatory state is the onset 
of vortex-shedding. Of the many examples, we cite here stall on airfoils when the angle of attack exceeds 
a critical value (Refs. 19»20) and the wake of the flow past a cylinder when the Reynolds number exceeds 
50 (Ref. 21). The equilibrium state resulting from a Hopf bifurcation Is often succeeded by bifurcation 
to a quasl-periodic equilibrium state with further increase of the parameter. To accommodate this and 
subsequent bifurcations, the integration regime in Eq. (22) must be interrupted at each such occurrence, 
and an Interval of 0(e) Inserted In which a discrete change, such as Eq. (23), to a new equilibrium state 
is allowed. The qua3i-periodic state may be succeeded by an aperiodic (chaotic) equilibrium state. As 
noted in sketch (e) In Fig. 4, appearance of the latter state signals the presence of a "strange attrac- 
tor." The slender body of revolution may undergo such a sequence of bifurcations as angle of attack is 
increased incrementally toward 90°. This example, among others, is discussed at length in Ref. 1. Sev- 
eral additional examples of bifurcation sequences relevant to aerodynamic applications are contained in 
Refs. 2 and 22. 

The last subdivision in Fig. 4, sketch (f), representative of aerodynamic systems having strong rate 
dependence, goes beyond the previous ones in that more is at Issue In the modeling than accommodating 
bifurcations of the equilibrium state. An adequate representation of the aerodynamic system In this cate- 
gory of problem represents a severe test for any mathematical model. We believe that the oscillating 
airfoil in what Is referred to as dynamic stall is exemplary In characterizing the nature of the problem 
in this category. Fortunately, a comprehensive set of experimental results exists for the problem, 
reported in a series of papers by McAlister, Carr, and McCroskey (Refs. 19,20). Included are excellent 
flow-visualization pictures which reveal many features of the physical mechanisms involved. In the 
following section we refer to these results and briefly demonstrate the extent to which our amended model, 
Eqs. (22) and (23), is capable of capturing their essentials. 

7.2 Dynamic Stall 

Figure 5, reproduced from Ref. 1, is representative of some of the principal results of Refs. 19 
and 20. Here, angle of attack a Is defined geometrically rather than kinematically as before. Flow- 
visualization pictures (obtained by means of the hydrogen-bubble technique, Ref. 19) of the airfoil oscil- 
lating about the quarter-chord clearly show that as the airfoil passes through a given angle of attack 
(a - 15°) with a increasing (a > 0), the flow has a very different character than it has as the airfoil 
passes through the same angle of attack with a decreasing (a < 0). The plots of C L versus a at the 
bottom of Fig. 5 (Ref. 20) reveal that the lift variations likewise are very different, depending on 
whether a is increasing or decreasing. With a increasing, the instantaneous value of the lift 



25-9 


coefficient can significantly exceed its static value (more precisely, in our terms, the mean value of the 

lift response in the equilibrium state), while the opposite is true with a decreasing. Finally, the 

differences between lift coefficient with a increasing and decreasing are dramatically enhanced as the 

reduced frequency k is increased incrementally. Increasing k may be interpreted^ increasing the 

magnitude of a, since with a - a + a sin u>t (Fig. 5) and k » WC/2V . we have |ac/2V I = a k. 

mo o 1 o 1 max o 

Thus, both sign and magnitude of a significantly influence the experimental results. These are 
important new features that a mathematical model ought to be capable of capturing. 

Let us first take up the issue of aerodynamic bifurcation in the equilibrium state for the airfoil 
under study. As the angle of attack i3 increased incrementally, a critical value of a is reached 
(a < 10°) at which separation of the flow begins, signaled by the appearance of a singular point in skin 
friction on the upper surface near the leading edge. With an infinitesimal increase in a, the singular 
point splits to form a half-saddle point of separation followed by a half-saddle point of reattachment, 
creating the familiar leading-edge separation bubble. The onset of separation is designated a structural 
bifurcation, inasmuch as it signals a change In the topological structure of the flow (Refs. 1,22). The 
change occurs uniquely (the original flow structure no longer represents a possible solution of the flow 
equations) and the flow remains steady. At an angle of attack of about 12°, an inflexional flow profile 
within the bubble becomes unstable, and periodic vortex-shedding begins. This is a Hopf bifurcation and 
available evidence suggests that it is supercritical. Figure 6 is a bifurcation diagram of the event with 
amplitude of the periodic lift fluctuation, designated C L (»;a(x)), as ordinate. The mean value of the 

lift coefficient in the equi librium sta te is shown in the same figure. To conform with the notation of 
Ref. 5, we designate it by C L («;a(x)). Beyond a - - 12°, the diversion of some of the available 

energy to maintain the fluctuation causes C L (»;a(x)) to begin to drop away from the values it would have 
taken (dashed curve) if the original stationary flow had remained stable. Its variation with a is that 
of the measurements for the curves labeled "static” in the plots of C L versus a in Fig. 5. In the 
absence of contrary evidence, we assume that no further aerodynamic bifurcations occur over the remaining 
range of angle of attack. Thus, Eqs. (22) and (23) should apply to the oscillatory motion illustrated in 
Fig. 5, with Eq. (23) required at a * a Q - 12°. 

Now let us examine the behavior of a typical lift response C L [a(£);t,x] for the part of a cycle in 
which a is increasing (a > 0), and let us choose a(x) to be somewhat larger than a =12°. During the 
portion of the motion (0 < £ < x) before a(£) is constrained to remain constant at aCx), the conditions 
correspond to those depicted via flow visualization in the left-hand photograph of Fig. 5. Let us focus 
attention on the flow component parallel to the airfoil surface at a station near the leading edge at 
time £ just prior to £ - x. We choose to observe the flow in the same reference frame as the one in 
which the camera photographing the flow wa3 situated. The origin of coordinates is fixed to the airfoil's 
axis of rotation, but, rather than rotating with the airfoil, the coordinates remain stationary relative 
to the oncoming wind. Within this reference frame, the instantaneous velocity profile at the fixed sta- 
tion may be viewed as having two components, as shown in Fig. 7(a). The first profile is the time- 
invariant mean profile which would exist alone were the airfoil to remain stationary at the instantaneous 
position a(£) = a(x). It has a region of reversed flow typical of the velocity profiles within the sepa- 
ration bubble in the equilibrium state. The second profile represents velocity induced by the motion of 
the airfoil. Its magnitude at the airfoil surface is directly proportional to the angular velocity a(£), 
and in particular, since £ * x, it is proportional to a(x). It is clear that with a(x) > 0, this con- 
tribution tends to diminish the region of separated flow, and more so as a(x) is increased. The situa- 
tion begins to change as soon as £ exceeds x, for then a Is constrained to remain constant at 
a(x). As shown in Fig. 7(b), while the first profile remains the same, the second profile immediately 
changes, since the velocity at the wall must go to zero as soon as the airfoil stops moving. With 
increasing values of £ this profile decays progressively, so that the zone of separated flow associated 
with the first profile gradually begins to appear. After a lapse of time, a separation profile emerges 
that is unstable to small perturbations. In response, the periodic fluctuation characteristic of vortex- 
shedding begins, and its amplitude grows toward the constant value corresponding to the amplitude of fluc- 
tuations in the equilibrium state at the chosen value of a(x) > a c . We show this schematically in 
Fig. 8. The interval of £ during which the onset of fluctuations i3 delayed is denoted by T; clearly, 
the duration of T is directly proportional to the magnitude of a(x). During the interval T, in which 
the flow over the upper surface of the airfoil remains more or less attached, the mean value of the lift 
coefficient tends toward a value representative of the attached-flow condition. A conservative estimate 
for it would be the value obtained from a computation based on inviscid flow theory; we indicate it on 
Fig. 8 as an upper bound (dashed line). After the fluctuation begins, the mean value of the lift coeffi- 
cient begins to decline toward the value it reaches in the equilibrium state. Thus, with a(x) > 0, the 
instantaneous mean value of the lift coefficient can overshoot its "static" value; the extent depends 
directly on the length of the interval T, which in turn depends directly on the magnitude of a(x). 

The behavior of the analogous lift response C^[a(£);t,x] for the part of a cycle in which a is 
decreasing (a < 0) is very much the reverse of the behavior just described. Let a(x) have the same value 

as before (a(x) > a c ), so that over the portion of the motion 0 < £ < x, conditions correspond to those 

depicted In the right-hand photograph of Fig. 5. In sharp contrast to the previous case with a > 0, the 
velocity profile induced by the airfoil's motion drastically extends the separated flow regime, particu- 
larly while the airfoil is moving, and the more so as a(x) is increased negatively. After £ exceeds 

x, whereupon a is constrained to remain constant at a(x), the violent fluctuations begin to die down as 

the second profile dies away. In contrast to the previous case, there is no interval in which the airfoil 
is free of fluctuations. The time it takes for the violent fluctuations to die down to the level 
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representative pf the equilibrium 3tate is directly proportional to the magnitude of |o(t)|. During this 
interval, the mean value of the lift coefficient must be^less than its "static” value; the extent of the 
difference again depends directly on the magnitude of |o(t)|. (We note in passing that large values of 
|a(t)| may represent a class where Fr£chet differentiability could depend on more than stability of flow 
in the equilibrium state. It may depend as well on the nature of the initial flow. The possibility will 
require further study.) 

This explanation of events suggested by the form of modeling we have adopted would appear to be in 
full accord with experimental observations. With one potential reservation, we conclude that a mathemati- 
cal model at the level of Eqs. (22) and (23) will be capable of incorporating the events we have 
described. 

7.3 Applications to Dynamical Systems 

We have seen that the amended mathematical model should be capable of describing the aerodynamic 
characteristics in all six forms of aerodynamic force (or moment) responses illustrated in Fig. 4. The 
forms involving aerodynamic bifurcations that result In time-dependent equilibrium states 
(sketches (d)-(f), Fig. 4) are of particular interest. Applied to problems involving dynamical maneuvers 
of aircraft, these forms will require extensive and novel experimentation. Problems in this category can 
be divided into two subclasses. 

The first subclass of problems involves oscillatory motions about an equilibrium state at which a 
bifurcation occurs. This was the case In the dynamic stall experiment studied In the last section.. Simi- 
lar cases need to be studied in which the oscillatory motions are free, rather than forced,, to determine 
how vortex-shedding frequencies may be modulated by the system’s natural frequencies, causing, for exam- 
ple, "frequency lock-in.” Another potentially important question is whether vortex-shedding (from, e.g., 
three-dimensional analogs of the airfoil’s leading-edge separation bubble) may be a source of wing rock.. 
Here, the phase relations between vortex-shedding on the two wing panels would appear to be a controlling 
factor. 

The second subclass of problems involves oscillatory motions that remain within a space free of crit- 
ical values of parameters. Here, simplifications in the modeling can be effected as previously (Ref. 8) 
in connection with time-invariant equilibrium states. Dynamical systems within this subclass of problems 
Involving vortex-shedding from obstacles have inspired a large collection of literature under the category 
of "vortex-induced oscillations.” Excellent surveys of the field have been published by Sarpkaya 
(Ref. 23) and very recently by Bearman (Ref. 24). The archetypal problem is the flexibly mounted cylinder 
Immersed in a uniform oncoming stream. Applied to this problem, our approach to modeling yields a differ- 
ential equation of the form 


h + 2ph + < 2 h - B(h)sin(w(h)t) (25) 

where the quantities h,y,tc characterizing the spring-mass system are illustrated in Fig. 9(a). The 
forcing term, which expresses the contribution of the fluctuating lift due to vortex-shedding, is novel in 
that the amplitude B and frequency o> are coupled to the velocity h on the left-hand" side. Our 
modeling approach provides an explanation as follows: Amplitude B is a function of instantaneous 

Reynolds number (Fig. 9(b)), which makes it a function of total velocity |/u 2 + h 2 . Strouhal number k, 
the dimensionless frequency of the fluctuating lift in the equilibrium state, is a constant, Independent 
of Reynolds number in the range of velocities of interest (Fig. 9(b)). This makes w,, the actual fre- 
quency, a function of h through 


k * » — * const. (26) 

tot 

We have found that Eq. (25) has a rich repertory of solutions, and appears to be capable of capturing 
distinctive features of the system’s behavior (frequency lock-in, amplitude Jumps, hysteresis) that care- 
ful experiments have revealed (Refs. 23,24). The equation’s richness is easier to understand if we trans- 
form it Into a set of autonomous first-order equations. There are three such equations, a number which 
suffices to permit solutions having "chaotic" behavior. We believe that study of simple systems such as 
this one can teach us how to perceive what otherwise might seem bizarre or even random behavior on the 
part of more complicated dynamical systems such as maneuvering aircraft. 


8. CONCLUDING REMARKS 

This paper is a continuation of our effort to provide a consistent formulation and theoretical method 
for studying nonlinear problems in flight dynamics. We have proposed a framework having parallel observa- 
tional and mathematical components to facilitate study of the aerodynamic contribution to the flight- 
dynamics system. The observational component consists of the elements structure, change, chaos, and 
scale. These are complemented by elements of a mathematical framework featuring the language of topology 
and concepts drawn from the theory of nonlinear dynamical systems. Here, particular emphasis has been 
placed on the element of change and its mathematical representation by means of bifurcation theory. 

Of the two principal realms of application of bifurcation theory to f light-dynamiC3 studies, we have 
concentrated attention on its role in the mathematical modeling of the aerodynamic contribution to the 
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aircraft's equations of motion. Our study stressed a number of themes which led to the following 
conclusions: 

1 . Studies in nonlinear flight dynamics can be enhanced by linking them to the common features of 
parallel studies in various allied fields. The unifying medium is the rapidly growing body of theory 
underlying research in nonlinear dynamical systems. Exemplifying this, our approach to modeling, involv- 
ing nonlinear indicial responses and generalized superposition integrals, was made compatible with one 
originating in the field of electrical circuits and systems, based on the use of nonlinear functional 
expansions. Operations involved in forming the indicial response were found to be equivalent to Fr^chet 
differentiation. Invalidation of Fr^chet differentiability proved to be the common element signaling the 
necessity of amending the mathematical model to accommodate potential discontinuous aerodynamic 
behavior^ Aerodynamic bifurcation is one means of invalidating Fr^chet differentiability. 

2. Aerodynamic bifurcation can be defined as the replacement of an unstable equilibrium flow by a 
new stable equilibrium flow at a critical value of a parameter. Acceptance of a list of presuppositions 
permitted extrapolating our understanding of aerodynamic bifurcations from a framework involving bounded 
laminar flows governed by the incompressible Navier-Stokes equations to one involving unbounded turbulent 
flows governed by modeled evolution equations based on the compressible Navier-Stokes equations. The 
issue of aerodynamic bifurcations represents a convergence of interests for the fields of turbulence 
modeling and the modeling of aerodynamic responses. Amended to accommodate aerodynamic bifurcations, with 
one potential reservation the mathematical model should be capable of describing the aerodynamic 
characteristics of all six major subdivisions that we have postulated for the aerodynamic response. 

3. Aerodynamic bifurcation i3 one of at least two means of invalidating Fr£chet differentiability of 
the aerodynamic response. A theory for enumerating these means offers a possible alternative to bifurca- 
tion theory. Inasmuch as it would incorporate the latter, it would be a potentially more inclusive way of 
classifying aerodynamic phenomena that are important in flight-dynamics applications. 
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Boundary conditions for lift at time t due to pulse in angle of attack at time C 



Fig. 2 Passage of angle of attack through a critical value a at which equilibrium flow becomes 

unstable. 









Fig. 6 Bifurcation diagram for amplitude of lift fluctuation. 



Fig. 7 Velocity profiles near leading edge, (a) £ < x. 


(b) £ > t» 
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Fig. 8 Mean and fluctuating lift responses. 




Fig. 9 Flexibly mounted cylinder, (a) Spring-mass system, (b) Bifurcation diagram for cylinder in 

uniform stream. 
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